Levitation of Current Carrying States in the Lattice Model for the Integer Quantum Hall Effect 
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The disorder driven quantum Hall to insulator transition is investigated for a two-dimensional lattice model. 
The Hall conductivity and the localization length are calculated numerically near the transition. For uncorrelated 
and weakly correlated disorder potentials the current carrying states are annihilated by the negative Chem states 
originating from the band center. In the presence of correlated disorder potentials with correlation length larger 
than approximately half the lattice constant the floating up of the critical states in energy without merging is 
observed. This behavior is similar to the levitation scenario proposed for the continuum model. 
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Soon after the discovery of the integer quantum Hall ef- 
fect (QHE) the fate of the current carrying states had been of 
much concern. At zero temperature and magnetic field B — 
a disordered two-dimensional electron gas should behave as 
an Anderson insulator. Therefore the question arose, how the 
transition from the QHE state to the insulating phase takes 
place when B goes to zero. Alternatively, for fixed magnetic 
field one can also consider the increase of the disorder and 
the broadening of the Landau bands. In the presence of dis- 
order, both Khmelnitzkii and Laughlin [|[] argued that the 
current carrying states that reside at the center of the Lan- 
dau bands cannot disappear discontinuously as _B — > 0, but 
rather float up in energy above the Fermi level. These gen- 
eral quasi-classical considerations which are based on scal- 
ing theory and supported by numerical results [^j] temporarily 
settled the case. Furthermore, the floating up of the critical 
energies to infinity became an essential cornerstone in con- 
structing a global phase diagram (GPD) [|| of the QHE. The 
levitation picture also implies that in the disorder -magnetic 
field plane only such transitions between QHE-states are pos- 
sible at which the change of the quantized Hall conductivity, 
A<Txy, is ±e^//i. Therefore, the final transition to the insula- 
tor must occur always from the lowest (i/ ~ 1) Hall plateau. 

Recently, the appealing levitation scenario has been ques- 
tioned both experimentally and theoretically and the debate 
about the fate of the current carrying states and the appearance 
of the phase diagram has been restarted [^-||]. One reason 
for this development is probably the fact that despite consid- 
erable efforts the microscopic origin of the levitation re- 
mains still unclear. Although transitions to the insulator from 
1/ = 1, as expected from the GPD, and the floating up of the 
current carrying states have been observed in some experi- 
ments [[sl p^pl] ], direct transitions from higher Hall plateaus 
(z/ > 2) which according to the GPD should not be possible 
have also been reported |^,|l2||T3],^ . On the numerical side, 
the observation of direct QHE to insulator transitions with 
AfJajj, > e^/hin tight-binding lattice models and the absence 



of floating states [|14|-|16|| have led these authors to conclude 
that due to lattice effects which are always present in real sys- 
tems the conventional floating up model which neglects the 
periodic background potential is inappropriate for explaining 



the experimental results. In addition, a new universality class 
has been assigned to direct transitions from the QH-liquid to 
the insulating phase where at the critical point the conductivi- 
ties are found to be equal, a^^ ~ a^^ [ |l2Hl6[ . 

In fact, the studies presented in showed that with 

increasing disorder the critical states do not float up in energy, 
but instead, the states with opposite Chern numbers from the 
center of the tight-binding band move downwards and finally 
annihilate the cuiTent caiTying states below the Fermi energy. 
In lattice systems the extended states can be characterized by a 
topological nonzero Chern integer [|l9|] . Even with next near- 
est neighbor hopping, which breaks the particle-hole symme- 
try, no floating up of states across the Landau gap, but only the 
moving outward of states with opposite Chern numbers, away 
from the center, has been found [^. This holds also for gen- 
eral magnetic fields commensurate with the lattice size, which 
are given by a rational number of flux quanta per plaquette, 
B ^ p/q ■ h/ (ea^) with p 1, that lead to Chern numbers 
for sub-bands which are ^ ±1 [^. Here p and q are coprime 
integers, and a is the lattice constant. 

Although successful in explaining the experimental obser- 
vation of direct transitions to the insulator from higher Hall 
plateaus, the lattice model investigated in possesses 
a serious drawback. If magnetic field and lattice effects are 
weak as in the experimental situation in GaAs heterostruc- 
tures an effective mass approximation should be appropriate 
and allow for the floating up results of the continuum model. 
In the latter, states with negative Chern numbers are absent so 
that the proposed annihilation cannot occur. 

The aim of our paper is to resolve this controversy and to 
reconcile the predictions for the continuum model with the re- 
sults of the lattice model. The key for achieving this goal is to 
refrain from unphysical realizations of the disorder potentials 
usually considered in this context. We will show that if the 
assumption of independent random site energies is abandoned 
and long range coiTelations of the disorder potentials are al- 
lowed as it is the case in the experiments [pljl, the floating 
up of critical states across the Landau gap is also seen in the 
lattice model. 

We consider a single-band tight-binding Hamiltonian de- 
scribing the properties of non-interacting particles in a disor- 



1 



dered two-dimensional system in the presence of a perpendic- 
ular magnetic field. 



H = 



En 



+ e 



), (1) 



<kl> 



where Ck is the fermionic operator on lattice site k. The 
magnetic field enters the transfer terms connecting nearest 
neighbors via the phases which in the Landau gauge read 
hki = ±27r(p/g)(ffc • ey)/a, if n = fk ± ae^, and bu = 
else, where is the position of site k and Cx, Sy are unit 
vectors pointing in the x and y directions. W is the disorder 
strength of the correlated random potentials Wk G [— M^, M^]. 
We choose V — 1 and p/q = 1/8 so that the tight-binding 
band splits into 8 sub-bands. 

The behaviour of the current carrying states is investigated 
through a numerical calculation of the localization length A m 



A^/(S,M^) = - lim i-hiTr|G+^p, 



(2) 



and the Hall conductivity 




axy{E,W;M) 



lim lim x 

e^0+ L^oc h LM 



(3) 



, L L 

7-1 ifi in ^ 



The are M x A/ sub-matrices of the advanced and re- 
tarded Green functions, G^ = {E — H ± ie)~^, acting in 
the subspace of the columns m, n on the lattice. The a;„ and 
yn are M x M matrices with elements {x„)ij — nSij and 
{yn)ij — i&ij- Periodic boundary conditions are applied in 
the y-direction for the calculation of \m while axy is obtained 
for Dirichlet boundary conditions. The length (x-direction) 
and width (j/-direction) are denoted by L and M, respectively, 
and e takes care of the proper thermodynamic limit. 

The task is accomplished by means of a recursive Green 
function method developed previously 23 1. This pro- 
cedure is well suited for calculating very long disordered 
2d systems of finite width without getting into trouble with 
storage-space problems. Therefore, also the correlated on- 
site disorder potentials Wk have to be generated iteratively 
site by site. This is implemented in analogy to the succes- 
sive creation of Ising spin chains. In this case the algo- 
rithm can be based on the Gibbs representation of Markov 
random fields [^]. The correlated random number Hxk,Vk 
at a given site k with coordinates Xk,yk is calculated with 
the help of two neighboring numbers already determined, 
^ixk.vk = {i^xk-a,yk +M:rfc,afc-a)/2 + 2 e"^ 7. The uncor- 
related random numbers 7 are drawn from an interval [—1,1] 
with constant probability density. If \Hxk,yk I > '^^e value is 
reflected back into the unit interval by fixk,yk ^ ±2 — i^xk^yk ■ 
The correlation strength is tuned via the correlation parameter 
G, where C ~ corresponds to the uncorrected case. The 
correlated random numbers ji generated in this way are mul- 
tiplied by the disorder strength W which results in the set of 
disorder potentials Wk G [—W, W] used in Eq. (jl]). 



FIG. 1. Grayscale plot of realizations for a 128 x 128 array of 
correlated disorder potentials with C = 0.0, 0.5, 1.0, and 1.5 



An example is shown in Fig. |l] where gray-scale plots 
for arrays (128 x 128) of correlated disorder potentials with 
G = 0.5, 1.0, and 1.5 are shown in comparison with the un- 
corrected case G = 0.0 (upper left). White spots correspond 
to disorder potentials Wk — W whereas black ones to Wk — 
—W. To better characterize the disorder potentials, we have 
calculated the correlation function iir(/9) = {wkWi) p^\f:^_p^\ 
averaged over all pairs of disorder potentials at sites which are 
a given distance \rk — ri \ apart. We find an exponential rela- 
tion, K{p) ^ exp{—p/r]{C)), where the spatial decay of the 
correlations is governed by the correlation length 77(G) which 
for example amounts to ?/(G) /a = 0.8, 1.1, and 1.9for corre- 
lation parameters, G = 0.5,0.7, and 1.0, respectively. Thus, 
for G < 0.7 the correlation length of the disorder potentials is 
less than the magnetic length. Is/a = 1.1284 for p/q = 1/8, 
while it is larger for G > 0.7. Our data follow nicely the re- 
lation ri{C) = 1/2 [ln(tanh(G))]~^ in analogy to Ising spin 
chains. 

First, to trace the QHE to insulator transition we calculate 
the Hall conductivity (Txy{W) at energy E/V — —1.5, which 
corresponds to a filling factor « 2, as a function of disor- 
der W. The result for various correlation parameters is shown 
in Fig. ^ Smooth transitions from axy = 2e^//i to zero are 
observed for correlation parameters G = 0, 0.5, 1.0 and 1.5. 
Similar direct transitions have been reported also for larger 
filling factors in the case of uncorrected disorder potentials 
|p^Jl7|,p||. This observation and the assignment of the crit- 



ical point to the disorder where axy = O.bve /h = Ux 



have led the authors of 



to conclude that for uncorrected 



2 



disorder potentials this transition corresponds to a new uni- 
versality class. Although the transition region gets broader 
with increasing correlations, no indication of a Hall plateau at 
(Jxy/ie^/h) = 1 can be observed at first sight. From this point 
of view our results of the Hall conductivity for correlated dis- 
order potentials seem qualitatively to be the same as for the 
uncorrelated case. 
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FIG. 2. The Hall conductivity versus disorder at _E = —1.5 (fill- 
ing factor ;/ ~ 2). The correlation parameter of the disorder poten- 
tials are C = (O), C = 0.5 (V), C = 1.0 (A), and C = 1.5 (O). 
The system length is L/a = 5 • lO** and the width is M/a = 128 
(C = 0.0), M/a = 160 (C = 0.5, 1.0), and M/a = 256 
(C= 1.5). 

Next, to find out the reason for the broadening of the tran- 
sition region in the correlated disorder case, we have calcu- 
lated the disorder dependence of the localization length. The 
position of the current carrying states can be extracted from 
the finite size scaling of the normalized localization length 
XMiE,W)/M. With increasing system width M, Xm/M 
goes to zero for localized states while it converges to a fi- 
nite value for critical states. The results for ly — 2 and 
C — 1.0 are shown in Fig. H where the averaged Am is dis- 
played versus disorder for different system widths in the range 
32 < M/a < 160. The averaging is performed over 3 to 5 
realizations of the correlated disorder potentials and, for the 
calculation of Aj\/, L is always larger than 5 • 10^ lattice spac- 
ings a. With increasing system width M two distinct peaks, 
which correspond to two separate divergences of the localiza- 



tion length at about Vl^i^' ~ 1.8 and Wc^' ~ 3.0, emerge 
out of the broad maximum. This implies that two transitions, 
from a„j = 2e^/h — > e^//i 0, must be hidden in the broad 
decay of the Hall conductivity shown in Fig. |[ 

In fact, a closer look at the Hall conductivity curves of 
Fig. ^ suggests that there is, hardly noticeable, a tiny shoul- 
der near a^y = e'^/h for correlation parameters C — 1.0 and 
C — 1.5. Taking the derivative daxy{W) / dW approximately 
by |Ao'j:j,(VF)/AV1^|, one sees that the peak corresponding 



.(2) 



to maximal steepness is not situated at 



xy — e /h, which 
should be the case for a direct transition, but instead near 
(Jxy/ie^/h) = 1.5, and a second somewhat weaker hump at 
about 0.5 can also be detected. This clearly means that there 



is no new single direct transition, but two distinct usual tran- 
sitions which, however, are not easy to discern in finite size 
studies. This difficulty of resolving the quantized Hall steps 
applies also to the experiments in addition to temperature and 
inelastic scattering effects which may mask the transitions. 
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FIG. 3. The normalized localization length versus disorder at fill- 
ing factor V = 2. The correlation parameter is C = 1.0 and the 
system sizes are M/a = 32 (□), M/a = 64 (O), M/a = 96 (O), 
and M/a = 160(A). 

To see the details of the QH to insulator transitions more 
clearly the energy and disorder dependence of \m {E, W) / M 
has been calculated in the filling factor range < < 2.5. 
Depending on the value of the correlation parameter we find 
the surprising result that with increasing disorder the current 
carrying states are either annihilated by the anti-Chern states 
or do float up in energy as seen in the continuum model. The 
essence of our extensive numerical calculations is represented 
in Fig. ^ which shows the location of the two lowest critical 
states within the disorder-energy plane which has been ex- 
tracted from the finite size scaling of Xm/M. 

We start with looking at the results for the correlation pa- 
rameter C = 0.2. With increasing disorder W the downward 
movement of the anti-Chern states (depicted by •) coming 
from the band center can be tracked. The last current carrying 
state, i.e. that from the lowest Landau band (marked by O) 
which moves slightly to smaller energies, disappears at a dis- 
order of about W/V « 2.5 when both Chern and anti-Chern 
states touch. This behavior is similar to the uncorrelated case 
(C ~ 0, Chern and anti-Chern states are labeled by □ and ■, 
respectively) where the last current cairying state vanishes at 
about W/V ~ 2.85, also shown in FigQ but hardly to discern 
from the C = 0.2 data in case of the Chern state. 

On the other hand, the floating of the energy Ec of the criti- 
cal states across the Landau levels without merging is seen for 
both C = 0.5 and C = 1.0. Here, the states move to higher 
energies when the disorder is increased. For a given energy 
that corresponds to a filling factor in the range 1.5 < < 2.5, 
the separation between the two critical disorders, which mark 
the current carrying states floating up in energy, is enlarged 
by enhancing the long range correlations of the disorder. The 
floating up is slower, i.e. a smaller shift in energy is caused 
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by the same increase in disorder, for larger correlation length. 
For C — 0.5, the states are closer and float up in energy within 
a small range of disorder values. These data are near a special 
correlation strength Cs{B) with rj ^ a/2 when the upward 
floating critical states and downward moving anti-Chern states 
meet halfway. In that situation the critical disorder Wc{E) is 
constant which implies that only an infinitesimal change in 
disorder is needed for the the lowest critical state to float up 
and for the anti-Chern state to float down. 
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FIG. 4. The position of the critical states in the energy-disorder 
plane for various correlation strength of the disorder potentials: 
C = 0.0 (□, the anti-Chern state is marked by ■), C = 0.2 (O, 
the anti-Chern state is •), C = 0.5 (O), and C = 1.0 (A). 

In conclusion, we have shown that floating up of current 
carrying states with increasing disorder can take place also in 
the lattice model, if correlation in the random disorder poten- 
tials is considered. In this case, no single direct transition from 
higher Hall plateaus to the insulator is seen. Instead, a one by 
one transition is observed as it is known from the continuum 
model. These usual transitions may be, however, difficult to 
detect both in experiments and in numerical finite size studies. 
We have shown that long range correlated disorder potentials 
will facilitate this task. Finally, by decreasing the correlation 
length, a crossover to the known annihilation scenario takes 
place where states with negative Chern number move down- 
wards and merge with positive Chern states. This crossover is 
found to occur at a correlation length of about half the lattice 
constant. 

This research is supported by the DFG-Schwerpunktpro- 
gramm "Quantum-HaU-Systeme" Grant No. SCHW795/1-1. 

Note added.-Aftei completing our work we became aware 
of a related paper [ p5[ ] where for Gaussian correlated disor- 
der potentials a floating and merging of extended levels was 



reported. The merging was not observed in our calculations. 
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